Phase transition in the fine structure constant 
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Within the context of mass- varying neutrinos, we construct a cosmological model that has a 
phase transition in the electromagnetic fine structure constant a at a redshift of 0.5. The model 
accommodates hints of a time variable a in quasar spectra and the nonobservance of such an effect at 
very low redshifts. ft is consistent with limits from the recombination and primordial nucleosynthesis 
eras and is free of instabilities. 



I. INTRODUCTION 

Measurements of the cosmic microwave background, 
large scale structure, the evolution of the Hubble param- 
eter from luminosity-redshift relation of type la super- 
novae along with the abundances of light elements in 
the universe strongly indicate the existence of a dark 
energy of unknown origin that acts against the pull of 
gravity The combined data favor an effective de- 
Sitter constant that nearly saturates the upper bound 
given by the present-day value (which we denote by a 
subscript to indicate redshift z = 0) of the Hubble 
parameter Hq « 10"'^^ eV. This yields a dark energy 
density: poE - "iMl^H^ - (2.4 x 10"^ eV)'', where 
Mpi ~ 2.4 X 10^^ GeV is the Planck mass. 

The coincidence of the neutrino mass scale with the 
dark energy mass scale is suggestive that there may be 
a link between these quantities. Measurements of atmo- 
spheric neutrinos have provided evidence (at > 15cr) for 

disappearing (likely converting to Vr) when propagat- 
ing over distances of order hundreds (or more) kilometers. 
The corresponding oscillation phase is consistent with 
being maximal and the oscillations require a neutrino 
mass-squared difference of Srn^^^-^ ~ 2.5 x 10~^ eV^ 2]. 
The disappearance oscillations have been confirmed 
by the KEK-to-Kamioka (K2K) and MINOS experiments 
over baselines of 250 km and 730 km, respectively. To 
implement a connection of dark energy and neutrino 
mass 9] in a concrete manner, Fardon, Nelson, and 
Weiner (FNW) H introduced a Yukawa coupling be- 
tween a sterile neutrino and a cosmic scalar field (dubbed 
the acceleron), such that the neutrino masses mjy. are 
generated by the vacuum expectation value A of this 
field, i.e., mv^{A). For simplicity hereafter we only con- 
sider a single nonvanishing neutrino mass, m^. The ac- 
tive neutrino mass is determined through a seesaw mech- 
anism by integrating out the heavy sterile neutrino with 
mass M{A), now correlated with the acceleron. This 
gives an effective potential 
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for regions in which nonrelativistic neutrinos dominate f4] 



for regions where relativistic neutrinos dominate. Here 
V is the fundamental acceleron potential, tud is a Dirac 
neutrino mass and {E^) is the average (relativistic) neu- 
trino energy. At the minimum of the potential, the neu- 
trino mass is determined in terms of its density Ui,. This 
creates neutrino mass dependence on the environment 
with possible relevance to solar and short-baseline 
neutrino oscillations as well as the cosmic neutrino 
background [3]. Interestingly, since the "constants" of 
the Standard Model depend nontrivially on the scalar 
neutrino density, my {A) could induce variations in the 
fine structure constant of quantum electrodynamics, a, 
which is the focus of this Letter. 

The fine structure constant has been measured in the 
spectra of distant quasars (QSO) for a number of ab- 
sorption systems. Early high redshift measurements of 
a with the Keck telescope found no discrepancy in com- 
parison with laboratory measurements of a to an accu- 
racy of a few parts in 10~* ^ and other observers also 
subsequently put upper limits on any discrepancy be- 
low the 10~^ level. However, a discrepancy of Aa/a = 
-0.57 ± 0.10 X lO"'^ was reported in Ref. Q. Further 
observations with the Very Large Telescope found no dis- 
crepancy at this level [l3|, but the parameter estimation 
methods are currently under debate pdj. Even if a dis- 
crepancy exists, it is not excluded that an effect may be 
imitated by a large change of isotope abundances over 
the last 10 billion years [l4|- Thus further observations 
are mandated to definitively decide whether or not a is 
truly constant. 

It has been proposed that a variation of a could result 
from the temporal evolution of a quintessence field ^1^. 
However, the model predicts a rather small variation of a 
from high redshifts to the present unless the quintessence 
field has unexpectedly undergone a rapid slowing in the 
recent past. 

Here we pursue the implications of mass-varying neu- 
trinos (MaVaNs) on the variation of a over cosmic time 
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scales. For a class of dependences of M{A) and V a tran- 
sition in the neutrino phase may occur as the neutrino 
density varies [13] ■ We show that this phase transition 
allows the existence of two distinct stable phases for a.^ 
Our paper is organized as follows. In Sec. |TT] we an- 
alyze the requirements that a MaVaN phase transition 
occurs and discuss the conditions under which the model 
can circumvent hydrodynamic instabilities that may be 
manifest in the nonrelativistic regime [3 . [l6j . Then, in 
Sec. Illli we study the corresponding implications for the 
time variation of the fine structure constant. We sum- 
marize in Sec. IIVI 



II. PHASE TRANSITION IN THE MASS 
VARIATION OF NEUTRINOS 

The stationary points for the potentials in Eqs. ([T]) and 
2) are given by 
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where V'{M) = dV{M)/dM in the two cases. We choose 
the phase of the singlet neutrino field so that M is real 
and nonnegative. 

To examine the possibilities for a multiphase structure, 
we note that in string-based discussions masses and cou- 
plings are determined by the minima of stabilized mod- 
uli [l7|. For our purposes, we take this to mean that 
M{A) ~ Mo[l + [A - Aoflf] in the vicinity of the 
stabilization point Ao^ where / is a positive constant. 
Without loss of generality, we set Ao = ^- With ^ 0, 
we make the assumption: (I) M{A) has a unique sta- 
tionary point at its absolute minimum Mq. (Although 
not essential for the discussion, we adopt the simplifying 
assumption that M is an even function of A.) As a con- 
sequence, both and have stationary points at 
^ = where dM/dA = 0. From Eqs. © and (g]) addi- 
tional stationary points will exist if the following permits 
solutions: 
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FIG. 1: Qualitative behavior of the effective potential. 
The dot-dashed line indicates the supercritical regime where 
riu > riv^, whereas the solid line indicates the subcritical 
regime where Ui, < Ui,^. 



To proceed, rather than examining the system in full gen- 
erality, we prefer to illustrate the possibilities by impos- 
ing a further condition, namely (II) M'^V'{M) is an in- 
creasing function of M . With the help of conditions I 
and II, Eq. ([5]) will have solutions in their separate do- 
mains if and only if rii, > riy^c- If this condition is ful- 
filled, the additional stationary points are the two mirror 
values ±,4min corresponding to the value M (assuming 
there is only one) for which there is a solution. These ar- 
guments then imply that the effective potential has the 
form in Fig. [U where the lower (upper) curve is valid for 

To illustrate these considerations, and to show how to 
circumvent the instability issues raised in Ref. [13] in the 
nonrelativistic reginie, we choose a slight variant of the 
original FNW form [J| for the acceleron potential. 



V[M{A)] = A" \n{\M{A)/Mo 



(7) 



where j = 2, 3 if i = NR, REL for the nonrelativistic and 
relativistic cases, respectively. Here 



normalized so that V = Q &i A ^ 0. Equation ([5]) be- 
comes 
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^ For a completely different mechanism that leads to an abrupt 
change in a, see Ref. llSll . 



where fc = 1, 2 for i = NR, REL, respectively. The critical 
neutrino densities may be calculated using Eq. In 
terms of the mass of the heaviest neutrino in a dilute 
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environment (toi/,o = itt-'d/Mo ^ \/ Sm'^^^^), we find 



(9) 



Now suppose that for nonrelativistic ly + the neutrino 
density is subcritical. Then the only stationary point is 
at M = Mo {A = 0) yielding a neutrino mass = rriijfi, 
which is independent of the neutrino density, so that 
there are no stability problems. Instability due to increas- 
ing densities can occur at a redshift at which both (a) 
n^, > n^_c and (b) the neutrinos are nonrelativistic, with 
temperature < 1.2 A [14]. However, recent work [l^ 
has shown that the instability can be avoided for suffi- 
ciently weak coupling of the neutrinos to the acceleron 
during the relevant cosmological era. We explicitly show 
below that there is indeed a window of late-time insta- 
bility for the FNW model Q , but that it can be avoided 
if the acceleron couplings are gravitational or stringy in 
origin. 

In the nonrelativistic regime, and for > n^^, 
from Eqs. ([S])-© we have = A'^/ui, with n^, = 
3C(3)/(27r^)r^. The neutrinos will be nonrelativistic if 
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< 1. 
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where (p^) is the mean square neutrino momentum in the 
Fermi-Dirac distribution. This yields T,y < 1.1 A, which 
effectively coincides with the criterion for instabilities in 



Ref. [lj|. With the condition 
instability is 



> n 
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Thus, instabilities may appear if A/m^^o ^ 0.23. 

Since = r^,o(l + z) (with T^^o - 1-7 x lO"'' eV), 
Eq. pT|) can be expressed in terms of redshift as 
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< 1 + z < 6.5 A_3 
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where A_3 = A/(10"3 gy). 

Instabilities may be avoided if the coupling between 
the acceleron and neutrinos (3 satisfies the inequality [l^ 



(3: 



din rUu 




d\\iM 


dA 




dA 



< 



CDM 



1 



2n„ 



Mpi 



(13) 



where ilcDM (S^u) is the cold dark matter (neutrino) 
dimensionless density. For values of interest, /3 < 
lO/Mpi So far our analysis has been independent 

of a particular choice of M{A); however, a test of this 
criterion, and in what follows, the results are somewhat 
dependent on this choice. We consider two simple forms 
(which also satisfy assumptions I and II above): 



M = MoC 



(14) 



and 



M = Mo cosh A/ f , 



(15) 



reminiscent of M-theory potentials for moduli 19]. We 
analyze the first case in detail, and only provide results 
for the second. From Eq. (|13p . 
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A bound on J^]// may be obtained by imposing the 
nonrelativistic criterion T^/A < 1.1. From Eqs. ^ and 



^2^j2 _ 3C(3) / 
° " 27r2 A I A 



which gives 
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T<V ^aT, 



(17) 
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In our discussion of the a variation we will find that 
A_3 ~ 0.6(TOiy,o/0.05 eV)^/^, so that for m^^o not much 
in excess of 0.05 eV, \A\/f < 1.7. From Eq. ^ we 
see that //Mpi > 0.34 serves as a sufficient condition to 
avoid instabilities. 

For the alternate M{A) of Eq. ([13, we find /3 = 
tanh(^//)l// < 1// for aU A, so that //Mpi > 0.1 
provides a sufficient condition for weak coupling and sta- 
bility. 



III. DISCONTINUITY IN THE FINE 
STRUCTURE CONSTANT 

Allowing for couplings between the acceleron and stan- 
dard model fields, 2 the free Lagrangian for the electro- 
magnetic field tensor i^^jy can be written as 



Am = - J ZpiA/Mpi) F^.F^"- 



(19) 



which on expansion about the present value of A, 
becomes 



Cam = (I+kAA/Mpi 



.)F^,F^^'' , (20) 



with = A — Ao and k = 9^^f|_4 • The field renor- 

1 /2 

malization A^^ — > A^/Zp to obtain a canonical kinetic 

1/2 

energy, generates an effective charge ejZp . Following 



^ We do not prescribe a mechanism (which may be desirable for 
technical naturalness) via which loop corrections to the acceleron 
potential are suppressed. 
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Ref. [13|, we expand to linear order about the present 
value Co, to obtain 



Aa 
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where a = e^/(47r). Equation (|2ip reflects our assump- 
tion that the variation in a is uniquely derived from the 
evolution of the acceleron. 

In order to accommodate the meteorite data pol |. 
which do not show evidence for a time-dependent a, we 
require that A not vary from ground state equilibrium 
{Ao = 0) for z < 0.5. Consequently, the model predicts 
no variation of a during this era, in agreement with ex- 
isting limits [2l| . 

From Eq. p2|) . for a transition at z = 0.5, 



A_ 



0.61(m^,o/0.05 eV)i/\ 



or equivalently 



PA 



4 X 10"^ 

PDE 0.05 eV 



(22) 



(23) 



This precludes A saturating the present dark energy. 
This ratio is roughly similar to the neutrino contribu- 
tion to the dark matter density. Large scale surveys and 
WMAP together constrain the neutrino energy density 
to be fi^ < 0.02, whereas terrestrial measurements of the 
neutrino mass indicate n„ > 7 x 10""*. 

From Eqs. ([T2|) and ([22]) . the density is supercritical 
and the neutrinos are nonrelativistic for 0.5 < z < 4. 
So there will be a tiny variation of a, in agreement with 
observations of absorption lines in the spectra of distant 
QSO [9]. The z-dependence of this variation may be 
obtained in a straightforward manner. From Eq. ([5]), at 
the minimum. 



M{A) 



l + z 

1 + 
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where Zc (— 0.5 in our case) is the redshift for which 
For the gaussian potential, we have 



'V.c ■ 



f 
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and for the cosh potential, 
M.eosh- 



l + z 

1 + Zc 



l + Zc 



(25) 



(26) 



At z = 2 (the intermediate point of the data), |.4|// ~ 
1.4(2.8) for the gaussian and cosh cases, respectively. 
Taking f/Mp\ > 0.34(0.1) for the two cases from the 
previous section, we find from Eq. (j21[) . 
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> 0.3 K 
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Thus, accommodating the possible variation of a at a 
level of 5 parts per million requires k ~ 10~^. 

The quantity k may be bounded by available limits 
on Aa/a during the eras of recombination (z ~ 1100) 
and big bang nucleosynthesis (BBN) (z ^ lO^*'). Since 
neutrinos are relativistic at these redshifts, there are no 
stability problems. It is straightforward to show that 
> n^^^ as soon as the neutrinos become relativistic. 
Then, the system is in the dot-dashed phase of Fig.[Tl and 
at the (mirror) minima the field A is given by relativistic 
case of Eq. ([S]), 
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yields 



Vln(lOz) 
cosh^^ (10 z) 



gaussian 
cosh 



(29) 



(30) 
(31) 



Inserting these in Eq. pT|) . we find for the recombination 
and BBN eras, 
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where b^cc = 3, 10 and 6bbn = 5, 26 for the gaussian 
and cosh potentials, respectively. (The slow variation 
with m^Q has been ignored in both cases.) Existing 
limits on the variation of the fine structure constant, 
|Aa/a| < 0.02 (at the 95% C. L.) for both recombi- 
nation [22] and BBN [2^, in conjunction with the lower 
bounds on //Mpi can now be translated into bounds on 
the coupling constant: 



K < 0.02 
K < 0.01 



recombination 
BBN , 



for the gaussian potential and 



K < 0.02 
K < 0.008 



recombination 
BBN , 



(33) 
(34) 



(35) 
(36) 



for the cosh potential. The bounds on k are a few orders 
of magnitude larger than the required value to accomo- 
date existing data in the subrelativistic regime. 



IV. SUMMARY 

In the MaVaN framework, we have constructed a cos- 
mological model that can accommodate limits on the 
variation of the fine structure constant on short time 
scales, as well as the potential observation of a varia- 
tion in a from distant QSOs. The model has a phase 
transition in the neutrino mass at a redshift z — 0.5 that 
gives a phase transition in a. The existence of this phase 
transition precludes that the vacuum energy associated 
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with the acceleron field saturates the present dark energy. 
To circumvent hydrodynamic instabilities we assumed a 
sufficiently weak coupling (perhaps stringy in origin) of 
the neutrinos to the acceleron during the cosmological 
evolution. The model is consistent with limits on Aa/a 
from recombination and primordial nucleosynthesis. 
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